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Abstract 

We prove an upper bound for the e-mixing time of the symmetric exclusion process 
on any graph G, with any feasible number of particles. Our estimate is proportional to 
Trw(g) lii(|V^|/e), where \V\ is the number of vertices in G and Trw(g) is the 1/4-mixing 
time of the corresponding single-particle random walk. This bound implies new results 
C^ . for symmetric exclusion on expanders, percolation clusters, the giant component of the 

Erdos-Renyi random graph and Poisson point processes in R''. Our technical tools 
include a variant of Morris' chameleon process. 



(N 

Q^ ■ 1 Introduction 

^^ ■ The symmetric exclusion process is a continuous-time Markov Chain defined on a weighted 

graph G = {y,E,{we}e<^E)-, where F is a set of vertices, ii^ is a set of edges and to each 
e £ E we assign a positive weight We > 0. For k < \V\, /c-particle symmetric exclusion on 



f^ ' G has the following informal description. 



Informal description ofEX{k, G): Start with k indistinguishable particles placed on distinct 
vertices of V. Each particle moves independently according to the symmetric transition 
rates given by the edge weights, except that moves to occupied sites are suppressed. 

This is one of the most basic and best studied processes in the literature on Interacting 
Particle Systems [151 116] . Literally hundreds of papers have been written on this process, 
but most of these results apply only to restricted classes of infinite graphs, such as the 
lattices Z'^. 

Exclusion processes over finite graphs have also been a testbed for the quantitative 
analysis of finite Markov chains. Coupling [1], comparison arguments [8], the martingale 
method for log-Sobolev inequalities [12^ I23j. and variants of the evolving sets technology 
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|201 [2T] have been variously applied to this process. Sharp results are known for some 
special cases, such as the complete graph [12] and discrete tori {'L/L'LY [20| [23]. 

In this paper we consider EX(A;, G) over an arbitrary finite graph and bound its mixing 
time in terms of the corresponding single-particle random walk, which we denote by RW(G). 
Our result is very general, but we will see that it nearly matches previously known mixing 
results for EX(A;, G) for very specific G and also gives new results in many interesting classes 
of examples. We will also argue that the kind of result presented here is of conceptual 
interest. 

1.1 The main result, and why it is interesting 

Recall that the e-mixing-time of an irreducible continuous-time Markov chain Q on a fi- 
nite set 5", with transition probabilities {qt{s,s')}s^s'£S,t>Oi and stationary (equilibrium) 
distribution vr, is given by the formula: 

TQ(e) = inf jt>0 : max(iTv(gt(s, •), ^) < 4 , (l-l-l) 

where cItv is the total-variation distance (cf. (|2.2.ip ). The 1/4 mixing time Tq(1/4) will 
also be called the mixing time of Q. Our main result is: 

Theorem 1.1 (Main result; proven in Section II. 4p There exists a universal constant 
C > such for all e G (0,1/2), all connected weighted graphs G = {V,E,{we}e£E) with 
\V\ > 2 and alike {1,...,\V\ -I}, 

TEX(fc,G)(e) < C ln(|l/|/6)TRW(G)(l/4). 

Our bound follows quite naturally if one assumes (heuristically) that the mixing time 
of EX(A:, G) is not much larger than that of k independent random walks on G, a process 
we denote by R\N{k, G) in what follows. 

[Heuristic assumption] T^x{k,G)i^) ^ C'oTRvv{fc,G)(£); C'o > universal. 

This assumption, if at all true, is well beyond the reach of present techniques. However, it 
is at least plausible, given that RW(fc, G) and EX(A;, G) are similar. 

It can be shown that TRvv(A:,G)(f) and TRvv(G)(f/^) are of the same order if e/k <C 1, 
thus our assumption is equivalent to: 

[Heurisitic assumption]' T^x{k,G){^) ^ C*! TRW(G)(f/^)) C*! > universal.. 

Recall the general inequality "Trw(g)('5) < C2 ln{l/6) Trvv(g)(1/4)", with C2 > universal, 
which is valid for any < S < 1/2 [T]. Applying this to our assumption, we obtain: 

[Heuristic conclusion] T^x(k.G)i^) ^ ^3 li^(^/f) '''rw(G)(1/4), ^3 > universal. 



Theorem 11.11 coincides with this for A;>|y|'^,c>Oa universal constant; whereas for other 
A; it is a strictly weaker result. 

We emphasize that what we just presented is not a rigorous proof of Theorem ll.il since 
we offer no good grounds for our heuristic assumption. What is interesting is that the 
Theorem does give an a posteriori justification for a weakened form of the assumption. We 
note that the bound "TEx(A;,G)(f) ^ C'Trvv(g)(1/4) ln(A:/e)" is tight up to constant factors 
for some G (e.g. discrete tori (Z/LZ) , d fixed |20j); therefore, in some sense Theorem ll.il 
is quite close to the best that one might hope for. 

Many other complex Markov chains are built from simpler processes that interact; ex- 
amples appear in eg. [DEI US]. Given our main result, it seems reasonable that, at least in 
some cases, the mixing time of these complex processes may be bounded in terms of their 
constituent parts. Some of the techniques we use to prove Theorem 11.11 are very specific 
to EX(A;, G), but it may be that some of the same ideas will turn out to be useful in other 
cases 

1.2 Connections with Aldous' conjecture 

Another motivation for our paper is a conjecture of Aldous' for the interchange process, 
which was recently proved in p]. The interchange process on G with k <\V\ particles can 
be informally described as follows: 

Informal description of \P{k, G): Start with k distinct vertices of V labelled 1,2, ... ,k all 
remaining vertices (if any) are labelled "empty". For each edge e, switch the labels of the 
endpoints of e at rate We- 

One can obtain EX{k,G) from \P{k,G) by "forgetting" the labels of the k particles. 
In particular, the contraction principle fT] implies that T^x{k,G)i^) ^ '^\P{k,G){^) fo^ ^^ 
1 < fc< |y| -1 and ah eG (0, 1). 

Aldous conjectured - and Caputo et al. recently proved [6] (see also [10]) - that IP(fc, G) 
and EX{k,G) always have the same spectral gap as RW(G) (or R\N{k,G)). This is a re- 
markable result, but it does not say much about the mixing times of these processes, since 
the bounds for T|p(;j(3)(e) or T^x{k,G)i^) that can be obtained from the spectral gap are 
typically very loose. 

Theorem 11.11 gives tighter relations between these mixing times. In the proof of the 
Theorem, we will show that the bound claimed for '^EX(k,G)i^) i^ the Theorem statement 
in fact holds for T|p(fc (j)(e) whenever k < \V\/2. Our proofs can be adapted to show that: 

Va e (0,1), 3Ca > 0, VG, VA: < a\V\ : J^p^k,G){e) < C7„ Trw(g)(1/4) H\V\/^)- 

That is, one can get a bound similar to Theorem 11.11 also for the interchange process, as 
long as the fraction of empty sites is bounded away from 0. Unfortunately, this leaves out 



Example 


Bound for TEx(fc,G)(l/4) 


{JLJULY with nearest-neighbor bonds pO] 

Typical largest percolation cluster in (Z/LZ)'^ [H [22] 

Typical Poisson process, in [0, Vf [5] (case a > d) 

Bounded-degree expanders 

Giant component of G„ ^/n- c > 1 |llj 


\Y\^l'^\u\V\ 
|y|2/^ln|y| 

|y|2/'^in|y| 

In^l^l 
In^lFl 



Table 1: Mixing time bounds for EX(A:, G) via Theorem ll.il in examples where no previous bound 
was explicit in the literature. We take d as a fixed parameter and assume fc ~ |V^|/2. 

the most interesting case of IP(|y|, G), which is a random walk by random transpositions in 
the group of permutations of V . Fortunately, the restriction on k does not make a difference 
for the exclusion process. 



1.3 Applications and comparison with previous results 

It is not hard to apply Theorem 11.11 to specific examples: all one needs is a bound for 
the mixing time of simple random walk on the given graph, and RW(G) is typically much 
easier to analyse than EX(/i;,G). The only example where we know Theorem 11.11 gives a 
suboptimal bound is in the case G = {Tj/LI) with the usual bonds, where the optimal 
bound, obtained by Morris [20] is of the order L^ In k whereas ours is about L^ In L (both 
for d fixed). Notice that this difference is only relevant for quite small k. 

Table [1] presents the bounds given by Theorem II. II in examples where no previous bound 
appears explictly in the literature. The references are to papers where the mixing times of 
the corresponding graphs are computed. We consider only k ~ |^|/2 and omit constant 
factors. 

In fairness, we note that a combination of canonical paths, log Sobolev constants and 
comparison arguments could in principle be applied to examples. This method is discussed 
in Section [A] in the appendix. However, we note that: 

• To the best of our knowledge, no good canonical paths bounds have been worked out 
for the examples in Table [U and it might be hard or impossible to do so; 

• Even if such bounds were obtained, there are natural lower bounds for how good they 
can be (cf. Section [A]), and Theorem 11.11 is at least as good as these lower bounds, 
up to the constants (it is actually better by a In |y| factor in the case of expanders). 



1.4 Key steps of the proof 

Our proof of Theorem 1 1 . 1 1 can be broken into two main steps. We first show that IP(2, G) 
always has a mixing time comparable to RW(G)o 

Lemma 1.1 For any weighted graph G, 

T|P(2,G)(l/4)< 20000 Trw(G) (1/4). 

The remainder of the proof is much harder and consists of bootstrapping this first 
Lemma to a larger number of particles. 

Lemma 1.2 (Proven in Section 16. Ij) There exists a universal constant K > such 
that for all connected weighted graphs G = {V,E,{we}eeE), all e £ (0)1/2) o-nd all k G 

{l,---,\V\/2}, 

T|P(fc,G)(e)<^T|p(2,G)(l/4)ln(|y|/e). 

Before we continue, we show how Theorem II . II easily follows from the two Lemmas. 
Proof: [Proof of Theorem II. Ij Combining Lemma 11.21 with Lemma II. II gives: 

T|P(fc,G)(e) < CTrw(G)(1/4) H\V\/e) if e G (0,1/2) and k < \V\/2 

where G = 20000 X. The contraction principle [1] implies: 

TEX(fc,G)(e) < T|p(fe,G)(e) < C^Trw(g)(1/4) H\V\/e) if e G (0,1/2) and k < \V\/2. 

However, EX{k,G) and EX(|y| — k,G) are the same process with the roles of empty and 
occupied sites reversed. In particular, TEx(A;.G)(f) = '^EX{\v\-k.G)i^) ^^r all e. □ 

We now give an overview of the main ideas involved in proving the two Lemmas. The 
proof of Lemma 11.11 relies on realizing that there are two classes of graphs. Some G are 
"easy" , in that two independent random walkers are likely to meet by time O (Tr\/\/(g) (1/4)) 
from any pair of initial states. In this case, an argument of Aldous and Fill [I] suffices to 
prove Lemma 1 1.11 (see Proposition I4.4|) . 

On the other hand, if G is not easy, then for most initial states two independent random 
walkers are very unlikely to meet by time (Trw(g) (1/4)) (cf. Proposition l4.5p . Intuitively, 
IP(2, G) and RW(2, G) are similar in the abscence of collisions, and we will use this to prove 
Lemma ll. II over non-easy graphs. The negative correlation property will be crucial for this 
part of the argument; see Remark 14.31 for details. 

The proof of Lemma 11.21 is considerably more involved. We first note that there are 
two methods in the literature for moving from mixing of pairs of particles to many more 
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particles, both of which were introduced by Morris |19tl20j. The first one [19] gives bounds 
for walks on the symmetric group by random transpositions. Unfortunately, the method 
seems to require too much from the process to be useful in our general setting. Moreover, 
the bounds given by that method would have a factor of ln(|y|) ln(l/e) where ours has a 
ln(|y|/e) term. 

Morris' other method was introduced in his study of symmetric exclusion over {1\UE)'^ 
[20j . The so-called chameleon process features particles that change color in a way that 
encodes the conditional distribution of the k-th. particle in IP(A;,G) given the other k — 1 
particles. It is this method that we will successfully adapt to prove Lemma 11.21 

One way to understand Morris' construction is that it reduces the analysis of mixing 
to the study of pairwise collisions between particles. The analysis for (ZyLZ)*^ is greatly 
facilitated by the explicit structure of the graph, something that we lack in general. This 
will require certain technical modifications of Morris' construction, of which we will try to 
make sense with remarks in our proofs. 

1.5 Organization 

Section [2] reviews some preliminary material. Section [3] discusses RW(G), EX{k,G) and 
\P{k,G), presents their joint graphical construction and reviews the negative correlation 
property. Section [4] presents the proof of Lemma 1 1.11 The chameleon process is introduced 
in Section [5l It is then used to prove Lemma [L2] in Section El but several lemmas are post- 
poned to Sections [7] to m It would be pointless to describe these steps now, but Section [6^ 
provides an outline of those sections. Finally, Section [10] presents some final remarks, and 
the Appendix contains some technical steps that are not particularly illuminating, 

2 Preliminaries 

2.1 Basic notation 

N = {0,1,2,3,...} is the set of non-negative integers and N+ = N\{0}. For n G N+, 
[n] = {i £ N4- : i < n} = {1, . . . , n}. If S* is a finite set, IS*! is the cardinality of S. For any 

ke[\S\], 

^^ {Acs : \A\ = k} 



k 
is the set of all size-/c subsets of S and 

{S)t = {s = (s(l), . . . , s{k)) e 5^= : Vi, j G [fc], "i ^ i" ^ "s(0 / s(j)"} 

is the set of all fc-tuples of distinct elements in S. 



Notational convention 2.1 The elements of {S)k will always he denoted by boldface let- 
ters such as X, with x(i) denoting the i-th coordinate of x. 



Notice tliat witli tliese symbols, 

S' 
k 



^f],\ish\ = i\s\h. 



A graph is a couple H = {V,E) where y / is the set of vertices and E C (2) is the 
set of edges. For each e £ E, the two elements a,b £ V such that e = {a, b} are called the 
endpoints of e. 

A weighted graph is a triple G = (V, E, {we}eeE) where {V, E) is a graph and We > 0, 
the weight of edge e, is positive for each e G E. When a graph G is introduced without 
explicitly defining the edge weights, we will assume that they are all equal to 1. We will 
assume throughout this paper that all graphs we consider are connected. 

2.2 Basic probabilistic concepts 

C[X] denotes the law or distribution of the random variable X. 

Given two probability distributions /U, z^ over the same finite set S, the total variation 
distance between them is given by several equivalent formulas: 

dTv{n,iy) = max{i^{A) - u{A)) (2.2.1) 

fdu (2.2.2) 

(2.2.3) 

(2.2.4) 

^ ses 

Another equivalent definition of d^y is 

dTv{n,u 

where the infimum is over all pairs {X,Y) of S*- valued random variables with C[X] = fi 
and C[Y] = u (such a pair is called a coupling of (/i, i^)). This implies that for any pair of 
S"- valued random variables X, Y defined over the same probability space: 

dTv{C[X],C[Y])<r{X^Y) 



ra.ax{ii(A) - 


- i^iA)) 


sup / „ 

/:5^[0,1] J 


fdfi- 


"^(Ks) - 


^(s))+ 


s€S 
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)=infP(X/y), 



We will need the following simple fact: if (for i = 1, 2) //j, f j are probability distributions 
on the finite set Si, 

dTvifJ-i X ^2, 1^1 X U2) < dTYifii,iyi) + dTv(/"2, ^^2)- (2.2.5) 

We will write Unif (S") for the uniform distribution on a set 5" 7^ 0. This is the normalized 
counting measure on S, if S is finite, or normalized Lebesgue measure over S*, if S* C M . 

2.3 Markov chains and mixing times 

For our purposes it is convenient to define a continous-time Markov chain over a finite set 
S" as a family of processes 

{{xnt>o ■■ seS} 

defined on the same probability space, with the following properties: 

1. For each s £ S, Xq = s almost surely; 

2. Each X^ is a "cadlag" path over S: there exists a divergent sequence 

To = < Tl < T2 < . . . 

and a sequence {sj}i>o C S with sq = s with X^ = Si over each interval [r^jTi+i). 

3. For each /i > and each cadlag path {xu)u>o taking values in S (in the sense of 2.), 

p (Xf_^^ = s' I Xf, = xt/, < t' < t) = P {XI' = s') almost surely. 

The last property is the so-called Markov property. It also implies that the law of 
(Xj*_^^)/j>o equals that of {X^')h->Q under the above conditioning. It is well-kwown that any 
such process is uniquely defined by its transition rates 

^(,, ,') ^ lim nXl=^ GS^s^ s') 

e\o e 

or equivalenty by its generator: 

We will usually make no distinction between a Markov chain and its generator in our 
notation. 

In this paper we will only work with irreducible chains, ie. chains for which for all A C S 
with A ^ 9, S\A 7^ 0, there exist a £ A, b £ S\A with q{a, b) > 0. It is well-known that 



such Markov chains have a unique stationary distribution vr, ie. a distribution such that 
if s* is picked according to vr independently from the {Xf)t>Q^s<^s^ then £[Xj**] = vr for all 
t > 0. Moreover, 

Vs G S, dTy{C[X^],7r) \ as t ^ +00. 

(The symbol "\" denotes monotone convergence.) The e-mixing time of Q is thus defined 
as in the Introduction: 

Jnie) = inf{t > : maxdTv(>C[Xf],7r) < e} (e G (0, 1)). 

We will often need two elementary facts about Markov chains and their mixing times. 

Proposition 2.1 ([13], equation (4.36) in page 55) LetQ be a Markov chain on finite 
state space S. Then for all < e < l/2|j 

TQ(6)<riog2(l/e)lTQ(l/4). 

Proposition 2.2 ([Ij, Lemma 7 in Chapter 4) LetQ be a Markov chain on finite state 
space S with symmetric transition rates. Then vr is uniform over S and moreover, for all 
<e <l/2 andt> 2JQ{e), 

for all s, s' £ S, with the same notation introduced above. 

We also make the following convenient notational convention. 

Notational convention 2.2 By definition, for any cadlag path {xt)t>o there exists a di- 
vergent sequence to = < ti < t2 < . . . with xt constant over [tj,tj+i) for each i > 0. For 
t > 0, we define xt_ to be the state of xt immediately prior to time t. That is, 

f xu_^ , if t = ti for some i> 1: 
y xt otherwise. 

Notice that xt_ = Xf^g for all 5 > sufficiently small. 

3 Random walks, exclusion and interchange processes 

In this section we formally define the main Markov chains in this paper: RW(G), EX{k,G) 
and IP(fc, G). We also present the standard graphical construction for the three processes 
at the same time, and then discuss the negative correlation property for EX{k,G). The 
material in this section is quite classical: Liggett's books [151 IB] are basic references, and 
the manuscript by Aldous and Fill [Ij contains some additional facts on IP(A;, G) as well as 
a presentation that is somewhat closer in style to ours. 



^The result in [T3] is for discrete-time chains, but the proof trivially extends to continuous time. 



3.1 Definitions 

The three processes we are defined in terms of the same weighted graph G = (F, E^ {we}e&E) 
with V finite (cf. Section [2]). We will be implicitly assuming that G is connected, in which 
case one can easily show that the chains defined below are irreducible. It will be useful to 
define the transpositions: 

{h \i X = a 
a \i X = h 
X otherwise. 

We also write /e(A) = {/^(a) : a £ A} and /^(x) = (/e(x(2)))ti for A G Q and x G {V)k 
(respectively). 

Simple random walk on G, denoted by RW(G), is the continuous-time Markov chain 
with state space V and transition rates: 

^(^'^) = |o otherwise. «-'-)^(^)2)- 

We will also consider the process R\N{k, G) that corresponds to k such random walks per- 
formed simultaneously and independently. Since the transition rates of these process are 
also symmetric, it follows that the stationary distribution of RW(A;, G) is Unit (y ) for all 

A;GN+. 

The k-particle symmetric exclusion process on G, denoted by EX(A;, G), is the continuous- 
time Markov chain with state space (^) and transition rates: 

The transition rates are again symmetric, and the stationary distribution of EX(A;, G) is 
Unif((r 

The k-particle interchange process on G, denoted by \P{k, G), has state space {V)k- The 
transition rates of IP(A;, G) are given by: 

This process also has symmetric transition rates, and its stationary distribution is Unit {{V)k) 
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3.2 The standard graphical construction 

We now present the standard graphical construction of these three processes. Graphical 
constructions are standard tools in the study of Interacting Particle Systems pLSJ ^^^ ^^^ 
usually attributed to Harris in the literature. The basic construction presented here will 
be elaborated upon later in the paper; see Section [5l For brevity, we omit all proofs in this 
subsection. 

Set W = X^ee-E'^e- We need a marked Poisson process, ie. a pair of independent 
ingredients given as follows: 

1. A Poisson process V = {ti < T2 < T3 < . . . } C [0, +00) with rate W. 

2. An i.i.d. sequence of -E- valued random variables ("markings") {enjneN, with 

Vn e N, P (e„ = e) = We/W. 

Let < t < s < +00 be given. We define a random permutation lu.s] ■ V ^ V 
associated with the time interval (t, s] as follows: if P n (t, s] = 0, /((.,] is the identity map 
on y. If on the other hand 

V n (t, s] = {tj : m < j < n} / 0, 

we set /(t s] = fe„ ° /e„_i o • • • o /g^; ie. /(^ .,] is the composition of each transposition f^^ 
corresponding to Tj G (t, s], and the transpositions are composed in the order they appear. 
We also set It = /(o,t] for t > and /(t^^] =identity map over V. 

Remark 3.1 Strictly speaking, we should worry about what happens ifVr]{t,s] is infinite, 
or (more generally) some finite interval (a, 6] in [0, +00) has infinite intersection with V. 
However, since the probability of any of this holding is 0, we will simply ignore these issues. 

Notice the following simple properties: 

Proposition 3.1 (Proof omitted) For all <t < s < r, I(t,r] = ^(s,r] ° ^(t,s]- 

Proposition 3.2 (Proof omitted) For all < t < s < +00, £[/(t,s]] = ^[^(7 si]- 

Proposition 3.3 (Proof omitted) Let < to < ti < ^2 < ••• < tj.. Then the maps 
I(ti_i,ti]! 1 ^ "^ ^ ^; fl'^e independent. 

Notational convention 3.1 We "lift" the random maps Ft^s] ^^ permutations of {^j^ and 
{y)k, which we also denote by I(t,s]' 

/(MlW = (/M(x(l)),/(M](x(2)),...,/(t,.](x(A;))) (x G (y)^) . 
For brevity, we will often write x^,^^,x^ instead of It{x), It{A), It{'x.) (resp.). 

11 



The key property of the graphical construction is: 
Proposition 3.4 (Proof omitted) Let to ^ 0. Then: 

1. For each x G V, the process {I(to,t+to]i^)}t>o is a realization of RW(G) with initial 
state X. 

2. For each A G (^), the process {I{to,t+to\{^)}t>a is a realization ofE.X{k, G) with initial 
state A. 

3. For each x £ iy)k, the process {I(to,t+to]i^)}t>o is a realization of\P(k, G) with initial 
state X. 

3.3 The negative correlation property 

EX{k, G) enjoys important negative correlation properties. In this paper we only need a 
very special result, which is contained in any of [15 ^ 114 1 [3]. 

Lemma 3.1 Given A G (A, let {Af}t>Q be a realization of EX{k,G) starting from A. 
Then for all u G (y)2 - ie. for all distinct u(l), u(2) gV -, we have: 

P ({u(l) G Al] n {u(2) G AD) < P (u(l) G Ai) P (u(2) G Af) . 

Using the construction in the previous section, we can write the above inequality as: 

P ({/r'(u(l)) G A} n {/r^(u(2)) G ^}) < P (/r^u(l)) G a) P (/r^(u(2)) G A) . 

The following is then immediate from Proposition! 



Corollary 3.1 (Proof omitted) Given u G {V)2, let {uj }t>o be a realization of IP(2, G) 
starting from, u. Then for all A C V : 

P ({u[(l) G A} n {uf (2) G ^}) < P (u[(l) G ^) P (uf (2) G A) . 

4 The dynamics of pairs of particles 

The goal of this section is to prove Lenima[TTTJ We fix a weighted graph G = {V, E, {we}eeE) 
for the remainder of the section (and of the paper). The definitions of RW(G), RW(A;,G), 
EX(k, G) and \P{k, G) are all relative to this graph. 
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4.1 Some facts on RW(2,G') and IP(2,G) 

Much of this section will involve comparisons between IP(2, G) and RW(2, G). The following 
notational convention will be useful. 

Notational convention 4.1 Given x G V"^, {x^ = (x^(l),x^(2)) : t > 0} denotes a 
realization o/RW(2,G) from initial state x. That is, the trajectories o/ x^(l),x^(2) are 
independent realizations of RW(G) with respective initial states x(l),x(2). 

We collect several simple facts about RW(2, G) and IP(2, G) that we will need later on. 
The first one is obvious eg. from the graphical construction. 

Proposition 4.1 (Proof omitted) For i = 1,2, £[x^(i)] = £[x[(i)]. 

The next Proposition is a direct consequence of (|2.2.5p . 

Proposition 4.2 (Proof omitted) The mixing times o/RW(2,G) satisfy Tr\/\/(2,g)(£) — 
TRW(G)(e/2)- 

Proposition 4.3 Let k e 'N be given. Then Trw(2,g)(2~^) < (^ + 1)Trw(g)(1/4)- 

Proof: Follows from the previous Proposition combined with Proposition 12.11 □ 

The next Lemma has the following meaning. Suppose t is so large that x|^ is close to 
equilibrium. In this case, E [</>(x|^)] is close to the uniform average of (j) over V'^, for all 
mappings < < 1. The Lemma shows that E [0(x^)l cannot be much larger than that 
average. This will require the negative correlation property (cf. Corollarv l3.ip . 

Lemma 4.1 Let (p : V^ ^ [0, 1]. Then: 

VeG(O,l/16),Vt>TRw(G)(e),VxG(y)2,E[0(xf)] <8Vi + 9 ^ ^. 

Proof: Define the "good set" of all a G F with nearly uniform probability: 



Good = <; a e y : max 

i=l,2 



(xf (0 = ") - ^ 



- \v\S 



We will show towards the end of the proof that: 

P (xf Good^) < 8\/i, (4.1.1) 

which (since < </> < 1) implies: 

E |c/>(xf)I(^)^\G„„d2(x{)] < 8Ve, (4.1.2) 
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On the other hand, notice that: 

E[0(xf)lGood2(xf)] = Yl P(xf = (a(l),a(2))) </.(a) 

aG(Good2)n(y)2 

< Yl ^P(^'(l) e {a(l)'a(2)}> xf(2) e {a(l),a(2)}) 0(a) 

aG(Good2)n(y)2 

(Cor. EU < j; P(xf(l) G {a(l),a(2)}) P(xf(2) G {a(l),a(2)}) </.(a) 

aG(Good2)n(y)2 

(Prop. SH) = J] P(xf (1) G {a(l),a(2)}) P(xf (2) G {a(l),a(2)}) ,^(a) 

aG(Good2)n(y)2 



(a(i) G Good) < 



E 



ae(Good2)n(y)2 



2 + 



VfV 



1^1 y 



(a) 



(Vi<l/4) < 9 Y 



aGy2 



1^ 



2 ■ 



Combining this with (|4.1.2p finishes the proof, except for (j4.1.1|) . To prove that, we let 
Bad = l/\Good. Notice that: 



-y/e |Bad| 
\V\ 



< 



El 



aeV 



M 



1 



(1)=«)-^ 



V\ 



+ 



1 



(xf (2) = a)-j^ 



as each a G Bad contributes at least -y/e/|y| to the sum. But the RHS equals: 
dTv(/:[xf (l)],Unif(y)) + dTv(/:[xf (2)],Unif(y)) < 2e 

since t > TRvv(G)(f)- We deduce: 
Ve|Bad| 



\V\ 



< 2e, or equivalently |Good| > (1 - 2y/^) \V\. 



Moreover, P (xf (i) = a) > (1 - 2^/e)\V\~'^ for all a G Good, hence: 

P (xf (i) G Good) > 1^ (1 - 2V-e) > (1 - 2V-ef > 1 - 4^^. 



Inequality (j4.1.1|) now follows from: 

P (xf Good^) < P (x[(l) Good) + P (xf (2) Good) 
(Propositions!]) = P (xf (1) Good) + P (xf (2) Good) < 8a/^. 



D 
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4.2 When collisions are nearly as fast as mixing 

Recalling Notational convention I4.H we define the first meeting time M(x) of RW(2, G) 



started from x G F^ as the smallest to > such that x^j(l) = x|^(2) (this is a.s. finite by 



ergodicity). We will also write 

M>t(x) = mf{ho > : xf+,Jl) = xf+,^(2)} 

for the time until the first meeting after t (this is a "time-shifted" meeting time). 
The following definition will be crucial for our analysis. 

Definition 1 We say that a weighted graph G is easy if: 

sup P(M(x) > 20000 Trw(g) (1/4)) < 1/8. 
xey2 

We note that all long enough paths and cycles are examples of easy graphs. Non-easy 
graphs include (Z/LZ)'^ for d > 2 fixed and L sufficiently large, as well as large expander 
graphs. The next Proposition proves Lemma [1.1 1 for all easy graphs via a coupling argument 
due to Aldous and Fill. 

Proposition 4.4 Lem.ma \1.1\ holds for all easy weighted graphs. 

Proof: [Sketch] Given G, Aldous and Fill [ll Chapter 14, Section 5] construct a coupling of 
IP(|y|,G) started from two different states u, v. Letting {uj^,vp}t>o denote the coupled 
trajectories, the following property holds: for each 1 < i < \V\, up(i),vp(i) behave as 
independent random walks up to their first meeting time, which we denote by Mj. After 
this time Mi, up(i) = vp(i), ie. the two processes move together. This implies: 

Vt > 0, dTviC[ui]XWi]) < P (uf ^ vf ) < J^P (uf (i) ^ vf (i)) < J]P(M, > t) . 

It is easy to adapt this to a coupling of IP(2, G) starting from given x,y G {V)2, so that, if 
{xj^,yp}t>o denotes the coupled trajectories, we have: 

Vt > 0, dTY{C[xi],C[yi]) < P (Ml >t)+¥{M2>t). 

Now both Ml and M2 are the meeting times of independent random walkers on G, which 
shows that: 

Vt > sup tfTv(^[x[],£[yf]) < 2 sup P(M(z) > t) . 

x,yG(V)2 zeV2 

For t = 20000 Trvv(g) (1/4) and G easy, the RHS is < 1/4. By convexity, this implies that: 

1 



sup dTv(/:[xf],Unif((y)2))< 
xe(y)2 

In other words, T|p(2,G)(l/4) < 20000 Trw(g) (1/4). □ 
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Remark 4.1 The argument of Aldous and Fill actually proves Theorem \1.1\ for all easy 
graphs; see fH Chapter I4, Section 5] for details. 

4.3 Long time to meet in non-easy graphs 

We now consider what happens when IP(2,G) is performed on a graph that is not easy. 
Our first goal is to show that independent random walkers take a relatively long time to 
meet from most initial states in V . 

Proposition 4.5 Assume G = (y, i?, {wejees) is not easy. Then: 

-^ Y, P(M(v) < 20Trw(g)(1/4)) < -^. 

Remark 4.2 In general we cannot guarantee that P (M(v) < 20TR\yv((^)(l/4)) is uniformly 
small over all v G {y)2- In particular, the probability of collision from adjacent v(l), v(2) 
might be much greater than the above bound. 

Proof: [Proof of the Proposition] Set T = Trvv(G)(1/4)- Since G is not easy, there exists 
some X G V'^ with: 

P(M(x) > 20000 T) > 1/8. (4.3.1) 

Consider some k £ N. Using the Markov property and the notation introduced in Sec- 
tion [321 one can write: 

P(M(x) > 40fcr) = E [l{A/(x)>40(fc-i)T}P (M>40(fc-i)T(x) > 40r I ^fo(k-l)T) ■ 

The conditional probability in the RHS equals P (M(y) > 40T) for y = x^/^_-|^sy, hence: 

P (M(x) > 40A;T) < ( sup P (M(y) > 40T) ) P (M(x) > 40(A; - l)r) 

VyeV2 J 

(...induction...) < ( sup P (M(y) > 40T) ) . 

Applying this to k = 500 and using the bound in (j4.3.ip gives the following with room to 
spare: 

sup P(M(y) > 40r) > 8-1/5°° > e^^/soo > ^ 
yey2 500 
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Fix some y ^ V'^ achieving this supremuni. Notice that M(y) > 40T holds if and only 
if yR{l) / yf (2) for all < t < AOT. If that is the case, y^or+/.(l) ^ y^T+hC^) ^^ ^11 
< h < 20T. Using the Markov property as before, we see that: 

^ < P(M(y) > 40r) < P(M>20T(y) > 20r) = ^ F (yfoT = v) P(M(v) > 20T) . 

VGV2 



Moreover, by ([2X2]) . 

1^1 



E ^ (y^V = v) F (M(v) > 20r) < j; IP(^Av)>20r) ^ ^^^(^[yi?^^]^ unif (y2)). 



vev^2 vey2 

Hgiicg' 

497 , ,^, R , ^^ ,,^,2nn V- F (M(v) > 20T) 

vev2 I I 

We finish by noting that, by Proposition 14.31 20T > Trw(2,g)(2~^^)) hence: 

dTv(/:[yi,T],Unif(y2)) < 2-19 < ^, 

and therefore 

^ F(M(v) >20r) 496 _ ^ 

^o WV ~ 500 ~ ~ 125' 

n 



4.4 If meeting takes a long time, IP(2,G') and RW(2,G') are similar 

We have just shown that the meeting is unlikely to be smaller than 20Trv\/(g)(1/4) from 
most initial states. We now show that IP (2, G) is similar to RW(2, G) until the first meeting 
time. 

Proposition 4.6 For any x G {V)2 and s > 0, 

dTv(/:[xf],/:[xi])<F(M(x)<s). 

We will only need the following simple corollary (proof omitted) in what follows. 
Corollary 4.1 For any x,y G (1^)2 o-i^-d s > 0, 

dTv{C[^i],C[yi]) < P(M(x) < s) + P(M(y) < s) + dTv{C[^f],C[y^]). 
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Proof: [Proof of Proposition I4.6| We present a coupling of {x.(}t>o and {x^joo such 
that the two processes agree up to M(x). The Proposition then follows from the coupling 
characterization of (iTv(") "Oi cf. Section \T2[ 

Our coupling is given by a continuous-times Markov chain on S = {V)2 x V^ with 
transition rates given by q{-, ••). The state space can be split into two parts, A = {(z,z) : 
z G (1^)2} and its complement A*^. 

• Transition rule 1: The transition rates from any pair (x, y) G A'^ to any other pair 
in S are the same as those of independent realizations of RW(2, G) and IP(2, G). 

• Transition rule 2: The transition rates from a pair (x,x) G A are determined as 
follows: 

— Transition rule 2.1: For each e ^ E with |en {x(l),x(2)}| = 1, 

9((x,x),(/e(x),/e(x))) = We] 

- Transition rule 2.2: If e G E satisfies e = {x(l),x(2)}, 



g((x,x),(/e(x),(x(l),x(l))))=u;. 
g((x,x),(x,(x(2),x(2)))) = w, 






— Transition rule 2.3: All other potential transitions have rate 0. 

Inspection of the marginals reveals that this indeed gives a coupling of {x^}t>o and {x(}t>o 
when started from an initial state (x,x) G A. Moreover, the two processes can only differ 
after a transition has occurred according to rule 2.2. The first time when this happens is 
precisely the first meeting time of {x|^}j>o. □ 

4.5 Proof of the mixing time bound for IP(2,G') 

We now use the tools developed above in order to prove Lemma ll.ll 

Proof: [Proof of Lemma II. 1| The case of easy graphs is covered by Proposition 14.41 so 
assume G = {V, E,{we}eeE) is not easy. Let x,y be given and T = Trvv(g)(1/4)- Notice 
that for all A C {V)2, if {xj }j>o, {yj^}t>o are defined over the same probability space, 

P (xioT G ^) - IP (yioT e A) = E [P (xioT e A \ x^qt) - ^ (yW e A \ y^or)] 

< E [dTv(F (xioT e • I 4ot) ,1P (yioT e • | y^oT))] • 
Maximizing over A yields: 

dTv(/:[xioT],^[yioT]) < E [dTvir (x^ e • I x^) ,f (yW e • I y^or))] • (4.5.1) 
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By the Markov property and Corollarv l4.lt 

^Tv(IP (xior e • I x^ = v) ,P (yioT ^ • | yioT = w)) = c^tv(^[v2ot],'C[w^ot]) 
< P (M(v) < 20T) + P (M(w) < 20T) + (iTv(>C[v^or]> >C[w|^or])- 

Proposition 14.31 implies the third term in the RHS is < 2^^^ for any v, ■w. Using this in 
conjunction with (j4.5.ip . we obtain: 

dTv(/:[xioT], ^[yioT]) < E [H4ot)] + E [Hyior)] + 2"'' (4.5.2) 

where (/)(z) = P(M(z) < 20T). Notice that < </) < 1. We may apply Lemma O and the 
fact that 20T > Trw(g)(2"^°) (cf. Proposition [21D to deduce: 

r / / NT 7 v^ P(M(v) < 2or) 
E [H4ot)] < 2-7 + 9 ^ ^ i^p -. (4.5.3) 

Applying the same reasoning to i?^(y20T) ^^'^ plugging the results into (|4.5.2p . we obtain: 

dMC[^ioT],C[yioT]) < 18 E "^^"^^r^p""""^ + 2-^ + 2-^^ (4.5.4) 

Finally, we use the fact that G is not easy, combined with Proposition 14. 5^ to deduce: 

1 8 

dTv(/:[xioT],^[yioT]) < Y^ + 2-^ + 2-6 < 1/4 (4.5.5) 

with room to spare. By convexity, 

dTv(/:[xioT],Unif((F)2)) < 1/4, (4.5.6) 

Since this holds for all x G {V)2, we have T|p(2,G)(l/4) < 40T, which implies Lemma ll. II 
for non-easy graphs. □ 

Remark 4.3 The first inequality in |^.5.3p follows from Lemma \4-i\ which is a conse- 
quence of the negative correlation property (cf. Lemma \3.1\ and Corollary \3.1\) . This is the 
first crucial use we make of negative correlation in this paper. 

5 The chameleon process 

In the previous section we determined the order of magnitude of the mixing time of IP(2, G). 
Going beyond two particles will require an important additional idea that is based on Morris' 
paper [20) . His idea is to introduce the so-called chameleon process to keep track of the 
conditional distribution of one particle in IP(A;, G). We will need a different process, which 
will nevertheless call by the same name. 
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5.1 A modified graphical construction 

We will need consider a variant of the construction of IP(A:,G) presented in Section [57 
Consider three independent ingredients: 

1. A Poisson process V = {ti < T2 < T3 < . . . } C [0, +00) with rate 2W. 

2. An i.i.d. sequence of -E-valued random variables {e„}„gN, with P(e„ = e) = We/W. 

3. An i.i.d. sequence of coin flips {cn}nm with P(c„ = 1) = P(c„ = 0) = 1/2. 

Recall the definition of /e from Section [3.11 and set fl = fe, fe =the identity function. 
We modify the definition of the maps /(^ ,,] from Section 13.21 as follows: if "P n (t, s] = 0, 
/(t s] is the identity map, as before. Otherwise, 

V n {t, s] = {t„ < Tn+l < ■■■ < Tm} 

and we set: 

iit,s] = feZ°---°fe::iof!:. 

The thinning property of the Poisson process implies that {t„ : c„ = 1} is a Poisson process 
with rate W. One can use this to show that: 

Proposition 5.1 (Proof omitted) The joint distribution of the maps I(t,s], < t < s < 

+c«, is the same as in Section \3. SX 



5.2 The chameleon process 

The chameleon process is built on top of the modified graphical construction. The definition 
of the process will depend on a parameter T > which we call the phase length, for reasons 
that will become clear later on. 

Given y G {y)k-i-, let 0(y) = {y(l), . . . ,y(/c — 1)} denote the set of vertices that 
"occupied" by the coordinates of y. The chameleon process will be a continuous-time, 
time-inhomogeneous Markov chain with state space: 

Ck{V) = {(z,i?,P,VF) : z G {V)k-i; the sets 0(z),i?,P,VF partition V}. (5.2.1) 

Notice that we do allow any of the R, P, W to be empty in the above definition. For a given 
(z, R, P, W) £ Ck{V), it will be convenient to refer to the vertices in the sets 0(z), R, P, W 
as black, red, pink and white (respectively). Notice that any vertex v £ V will belong to 
one of these color classes. 

The evolution of the process from initial state (z, R, P, W) will be denoted by 

{izf,Rf,Pf,Wf)}t>o. 
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Box 5.1 The three kinds of updates in the chameleon process. 

• Constant-color phases: If t = t„ G {{2k - 2)T, {2k - 1)T] for some k £ N+, update: 

{zf,R?,pf,wf) = {re:{zi),r,:{Ri),re:{Pi),r,:{wf_)). (5.2.2) 

That is, the states of the endpoints of e„ are flipped if c^ = 1, and nothing happens 
if c„ = 0. 

• Color-changing phases. If i = r„ G {{2k — l)T,2kT], for k G N+, update as above 
unless: 

1- Sn = {w,r} has a white endpoint w G W[^ and a red endpoint r G Rf_; 
2. \P[:\<mm{\Rl\,\Wf_\}. 

If 1. and 2. hold, r and w both become pink and we call t a pinkening time. 

{z^,R^,Pf,Wf) = {zl,Rl\{r},Pf_U{r,w},Wf_\{w}). (5.2.3) 

• Depinking times. If t = 2kT with A; G N (t = or t lies at the end of a color-changing 
phase) and \P[^ \ > ram{\W[^ |, l-R^l |} (more pink than either white or red), flip a fair 
coin dk and make all pink particles become red or white depending on whether dk 
comes out heads or tails (resp.): 

{z,,R,,Pt ,W, )-| (^c^^c^0^p^CupC)^ ^^^0, (5.2.4) 



By definition, this process will only be updated at the times r„ (n G N) given by the Poisson 
process and at deterministic times 2kT, /c G N. Moreover, the updates at times Tn are of 
different kinds depending on whether r„ G {{2k — 2)T, {2k — 1)T] for some k G N+, or 
Tn G {{2k — 2)r, {2k — 1)T] for some k G N+. The three update rules are described in Box 
531 



Remark 5.1 Technically, this process is not cddldg, as it changes at time 0. However, 
one can see that it is cddldg outside of this point, so we will continue to use our notational 
convention 

{zi,Rl,p[:,wfj 

for all t > (cf. Notational convention \2.2\) . For t = 0, this simply denotes {z,R,P,W). 
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Remark 5.2 We briefly note that our chamaleon process is more complicated than that 
of Morris 120^ . In brief: his process does not have constant- color phases, and will depink 
right when the number of pink particles exceeds the mininiuTn of red and white. The second 
difference is a matter of convenience, but the first one will be fundamental at key steps of 
our argument. 

5.3 Two basic properties 

The next two results will be useful later on. We only sketch the proofs. 

Lemma 5.1 

{io,R^,Po,Wo) = {z,R,P,W) 

denote the initial state of the process and 

{zk,Rk, Pk, Wk) = the value of (zg.^_ , Rg^. , PScT. > ^2^T_ ) (^ ^ N+). 

Then {{z^, Rk, Pk,Wk)}kef>i is a discrete-time, time-homogeneous Markov chain. Moreover, 
if Dj is the j-th depinking time of the process, then Dj = Dj/2T is a stopping-time for this 
discrete-time Markov chain. 

Proof: (Sketch) Markovianity and time-homogeneity are obvious. To prove the stopping 
time property, it suffices to check that (setting Dq = 0), 

Vi>0, ^ = inf|/c>^ : \Pk\>mm{\Rk\,\Wk\} 
where we allow the inf to be +oo if the set is empty or -Dj-i = +oo. □ 

Lemma 5.2 Suppose (z^^^, i?§y,-P^'r, VF^^,) is the state of the chameleon process at time 
2iT (ie. at the beginning of a constant- color phase). Then: 

(^{2i+l)T' -^{2i+l)T' P{2i+l)T-'^{2i+l)T) = {H'^2iT) ^ HP'2iT) ^ HP2iT) ^ ^ 0^2iT)) 

where I = I(2iT,{2i+i)T] is the map defined in the modified graphical construction. 
Proof: By inspection. D 
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5.4 The chameleon process and conditional distributions 

We now explain the relationship between the chameleon process and conditional distribu- 
tions. 

Notational convention 5.1 x = (x(l), . . . ,x(A;)) G iV)k is represented as a pair (z,x), 
where z = (x(l), . . . ,x(A: — 1)) G {V)k-i and x = x(/c) € y\0(z). (Notice that Xj = {zj-,xj-) 
for all t > 0.) 

Proposition 5.2 (Proof omitted) Given an initial state x = (z,x) E (V')a; for \P{k,G), 
set R = {x}, P = $ and W = y\(0(z) U {x}). Consider the interchange process {x^ = 
(z|^,Xj)}j>o started from state x and the chameleon process {{z.^ ,R^ ,P^ ,Wi)}t>o started 
from configuration {z,R,P,W) G Ck{V). Then: 

yt > 0, Vb = (c, b) £ {V)k, P (xf = b) = E [inkt(6)I{^c=^|] (5.4.1) 

where 

ink,(t;) = I^,,^c} + i^ {v G V). (5.4.2) 

This is almost identical (up to changes in notation) to [20^ Lemma 1], and we omit its 
proof. It will be useful to think of '\nkt{v) as the amount of "red ink" at vertex v € V: a 
red vertex has one unit of red ink, a pink vertex has half a unit, and black or white vertices 
have no ink. We will see below that the total amount of red ink in the system determines 
the rate of convergence to equilibrium of \P{k,G). 

6 From 2 to A; particles via the chameleon process 

In this section we present the proof of Lemma 11.21 modulo several Lemmas about the 
chameleon process that we will prove later. We then outline the remainder of the paper. 

6.1 Proof of Lemma 11.21 

Proof: We assume we have defined a chameleon process over CkiV) as in Section [221 We 
will take the notation and definitions from that section for granted. We also define: 

\pC\ 
\nkt = ^\nkt{v) = \R?\ + ^-^. (t>0) (6.1.1) 

vev 
We note for later reference that 

inki= Y, '"'^tW (6-1-2) 

v€V\Oizl) 
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since the vertices in 0(z() have zero red ink. 

We have argued in Proposition 15.21 that the distribution of \P{k,G) started from x = 
(z, x) S {V)k corresponds to a chameleon process started from (z, {x},0, y\(0(z) U {x})). 
Letting inkf denote the value of inkj in that chameleon process, we will show that: 



Lemma 6.1 (Proven in Section 18.21) The following inequality holds for all 1 < k < 

\V\-1: 



sup dTviC[xl],Vmi{{V)k)) <2k sup E 



xG(V)fc 



inkf 



where 



Fil 



xG{V)fc 



lim inkf = \V\ - k + 1 



\V\-k + l 



Fil 



The main goal is to bound the expected value in the RHS of the inequality in Lemma f6.ll 
Fix some x G {V)k and let Dj{x) denote the j-th. depinking time for the chameleon process 
corresponding to X. Also set ink- = ink^/^-, for this process. We will show in Proposition l7.ll 
that there are infinitely many depinking times, ie. there are infinitely many times of the 
form 2kT at which the number of pink particles is at least as large as the minimum of the 
numbers of white and red. The definition of the chameleon process implies that inkf can 
only change at depinking times, hence for any i>Oinkf = lift< Di(x) and inkf = ink- 
if -Dj(x) < t < L'j+i(x) for some j. We deduce that: 



ink? 



ink^ 



^" m-fc+i -Si' \v\-k+i 



+1 



{D,(x)>t} 



< 



E^ 



ink^ 



+1 



m>j 



\V\-k + l 
Taking expectations, we see that the RHS of the inequality in Lemma |6. II is at most: 

-X 

ink^„ 



{D,{x)>t}- 



2k sup {J2^ 



1 



Fill 



+ F{Dj{x) > t I Fil 



.1.3) 



\V\-k + l 
A simple (but technical) proposition proven in the Appendix will take care of the first term. 



Proposition 6.1 (Proven in Section [B]) For all £ > 1 and x G iV)k, 



E 



ink« 



\V\-k + l 



Fil 



<V\v\-k + i{^^\ 
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We thus have: 



2k sup < 

xe(y)fc 



m>j 



ink^ 



V -k + 1 



Fill 



+ P(Dj(x) > t I Fil 



< C2 \V\^/^ e-"'^ + 10k sup P {Dj{x) > t \ Fill) , (6.1.4) 

xe(y)fc 



where ci = ln(72/71) > and C2 = 720 are universal constants. 

Boundmg P(Dj(x) > t \ Fill) is the key step in the proof. Up to now all of our results 
have been valid for all values of k, \V\ and of the phase length parameter T > 0. The next 
Lemma will require restrictions on these values. 



Lemma 6.2 (Proven in Section IQ.Sft There exist universal constants C^^d^ > such 
that, if \V\ > 300, T > C3 T|p(2,G)(l/4) and k/\V\ < 1/2, then: 



Vx G {V)k, Vj G N : E 



e c'42' iFil 



<e^ 



e-' . 



If |y| > 300 Markov's inequality allows one to deduce that, for yet another universal 
constant L = C3C4 

P(L»j(x) > t|Fill) < e^~^^iP(2.G)(i/4). 
Plugging this into (j6.1.4|) and Lemma |6. 11 we obtain: 

dTY{C[^i]),\Jmi{{V)k)) < Ci \Vf/^ e~''^ + 10|y| e^"^^iP(2,G)(V4) _ 
Since this inequality holds for all j, we can take 

t 



J 



.2i^T|p(2,G)(l/4) 



and obtain: 



dTy{C[^i],Vn\i{{V)k)) < Ko\Vf/^e 



3/2„ 2LT| 



IP(2,G) 



(1/4) 



with Kq > universal. Comparing with the definition of mixing time in (ll.l.lh and noting 
that Unif((y)fc) is stationary for \P{k,G) finishes the proof in the case \V\ > 300. 

The case \V\ < 300 - ie. \V\ bounded by a universal constant - can be dealt with in 
several ways. For example, one may use the result of Caputo et al. [6] for the spectral 
gap of \P{k,G) together with the standard lower bound for Trw(g)(1/4) in terms of the 
spectral gap and the usual upper bound for T|p(^, £")(€) in terms of its spectral gap (see eg. 
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|17j for these standard bounds). Alternatively, one may use the analysis of Aldous and Fill 
(see Remark I4.ip together with the inequality 



P(M(x)>2fcTRW(G)(l/4)) < 1 



1 -' 



A\V\ 
which one can prove via Proposition 12.21 and a few simple calculations. □ 

6.2 Outline of the missing steps 

We now summarize the main steps left in the proof. 

1. In Section [7] we collect several facts about the quantity ink. The proof of Proposi- 
tion [6rT] on the decay of E 1 — ink^/(|y| — A; + 1) | Fill , presented in the Appendix 
(see Section |B|), relies on results from this section. 

2. Section [8] contains the proof of Lemma l6.lt which is based on an auxiliary result on 
conditional distributions (Lemma [2D 



3. Section [9] bounds the right tail of the first depinking time in a chameleon process, and 
then uses this to bound the exponential moment of the j-th. depinking time. This 
leads to the key Lemma 16.21 proven in Section [9.31 

7 A miscellany of facts on ink 

In this section we prove several facts we will need about the quantity inkj introduced in 
(j6.1.ip . We will use the same notation introduced in the proof of Lemma 11.21 (cf. Sec- 
tion EH): 



1. x G {V)k is some fixed state; 

2. {z,R,P,W) = {z,{x},i/},V\{0{z) U {x})) E Ck{V) is the initial state corresponding 
to X in the sense of Proposition 15. 2| 



3. inkf is the total amount of ink in (zp, i?p, Pf, Wf) (with the above initial state); 

4. Dj{x) is the j'-th depinking time for this process; 

5. finally, ink- = ink^ 



-X 



We will mostly omit x from the notation in what follows. Most proofs in this section follow 
by inspection, so we will be quite brief. 

In principle the total number of number of depinking times could be finite. We begin 
by showing that this is not the case. 
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Proposition 7.1 The number of depinking times is alm,ost surely infinite. 

Remark 7.1 Notice that this only works because our definition of a depinking time allows 
for "trivial depinking times" where there are only red and black (or only white and black) 
particles left. This was noted in Box \5.1\ 



Proof: [Proof of Proposition 17. 1| We use the following simple fact (proof omitted) : there 
exists some 5 > ^ such that each color-changing phase that starts with min{|i?p|, |VF/^|} > 
will have a pinkening with probability > 5, regardless of the past. This implies that, given 
s > 0, the values of |P/^| for t G [s, +oo) will have a strictly positive probability of increasing 
by the end of each color-changing phase, at least until \Pt^\ > min{|i?p|, jVl^/^l}. Since \Pf\ 
can only decrease at depinking steps, this shows that \Pf\ must continue to increase until 
\Pf\ > min{|i?p|, |VF^*^|}, and the next time of the 2kT will be a depinking time. □ 

The next result follows by inspection. 

Proposition 7.2 (Proof omitted) < inkj < |F| — /c + 1 for all j G N, a.s.. 
We now compute the amount of change of ink in each step. 

Proposition 7.3 For j € N, inkj+i G {\nkj + A(inkj), inkj — A(inkj)} a.s., where: 

min{r, \V\ — k + 1 — r} 



Ar) 



(r G N). (7.0.1) 



Moreover, conditionally on {ink^};^^Q, each possibility is equally likely. 

Proof: Box 15.11 shows that there are no pink particles left in the system after depinking 
is performed. This implies that inkj = ink^). = |-R§.|. Moreover, since the total number 

of non-black particles is \V\ — k + 1, there must be inkj red and \V\ — k + 1 — \r\kj white 
particles at time Dj. 

A pinkening step decreases the number of red and white particles by 1 each and increases 
the number of pink particles by 2. However, no pinkenings are performed if the number 
of pink particles is at least the number of red or the number of white particles. In other 
words, the number of pinkening steps until the next depinking is precisely the smallest p 
satisfying 2p > inkj — p or 2p > \V\ — k — 1 — inkj — p, which is p = A(inkj) for A defined 
in (|7lU]l . 

At the depinking step, the pink particles either all become white, or they all become 
red. These possibilities corresponds to inkj+i = inkj — A(inkj) or inkj+i = inkj + A(inkj), 
respectively. Which possibility will actually occur depends on the value of the fair coin 
dk that is flipped at the depinking time 2kT = -Dj+i- it is easy to see that the coin is 
independent of {\nki}i<j, and this implies that both possibilities are equally likely. □ 
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The next leinnia summarizes the above sequence of propositions and adds a useful 
remark. 

Lemma 7.1 The sequence {inkj}j>o is a Markov chain with initial state inkg = 1, absorb- 
ing states at and \V\ — k + 1 and transition probabilities given by: 

p{a, b) = - (I{6=a+A(a)} + I{fe=a-A(a)}) (a, b G {0,1,2, ... ,\V\ - k + 1}). (7.0.2) 

Moreover, it is almost surely absorbed in finite time in either or \V\ — k + 1. Finally, the 
event: 



Fill 



has probability 1/{\V\ — k + 1). 



hm \nkj = \V\ - k + 1 



j->-+oo 



(7.0.3) 



Remark 7.2 The event Fill corresponds to the number of red particles converging to \V\ — 
k + 1, i.e. that there are only black and red particles at all large enough times, or, equiva- 
lently, to red ink filling up all available space. Notice that we can rewrite: 



Fill = <! Hm inki = \V\-k + l\ , 

t— >-+oo 



which is the form that appears in the proof of Lemma \1.1[ 

Proof: [Proof of Lemma[7?T] The first sentence is obvious given the sequence of Propositions; 
only notice that p{a,a) = 1 if a G {0, \V\ — k + 1}. We omit the trivial proof of the next 
assertion, which implies inkoo = limj_j._|_oo inkj G {0, \V\ — k + 1}. 

Now notice that the increments of inkj are unbiased; that implies that this process is 
also a martingale. We thus have: 



n 



P(Fill) 



(inkooe{0,|F|-A; + l}) 



({inkjIjgN bounded, cf. Prop. 17. 2p 



(mart, property + in kg = 1) 



inkr. 



\V\ -k + l 



E 



inkf- 



\V\-k + l 



limj^+oo IE 



ink. 



\V\-k + l 



E 



inkfi 



\V\-k + l \V\-k + l' 
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We will need one final lemma before we proceed. 
Lemma 7.2 For all b G {V)k-i and t > 0, 

P ({zp = b} n Fill) 



,c 



b) 



\V\-k + l 



Proof: This follows from the previous Lemma if we can show that Fill and zp are inde- 
pendent. To see this, simply notice that Fill is entirely determined by the coin flips dk 
performed at the depinking times, whereas the value of zp does not at all depend on these 
coin flips. D 

Remark 7.3 It transpires from the above that the chameleon process conditioned on Fill is 
the same as the unconditional process, except that the coin flips dk performed at depinking 
times are biased. This remark will be useful in the proof of Lemma \6.S\ in Section \9.3X 



8 Convergence to stationarity in terms of ink 

In this section we will prove Lemma 16. H used in the proof of Lemma 11.2! (cf . Section 16. ip , 
in which we show that the amount of ink in the system can be used to bound the distance 
to the stationary distribution. We start with a preliminary result on marginals. 

8.1 The convergence to equilibrium of conditional distributions 

We will again use Notational convention 15.11 whereby any x = (x(l), . . . ,x(/i;)) € {y)k is 
written as a pair x = (z,x) with z = (x(l), . . . ,x(A: — 1)) and x = x{k). 

Let X = {z,x) G {V)k and consider the IP(A;,G) process {xj}t>o. Set R = {x}, P = % 
and W = y\(0(z) U {x}) and recall from Proposition 15.21 that the chameleon process 
{(zf ,i?f,Pf ,Wf )}t>o satisfies: 



Vt > 0, Vb = (c, b) G (y)fc, P (xf = b) 



E 



i{zf=c} 



inkf(6) 



(8.L1) 



where (as before) we use inkf (•) to denote the amount of ink in this chameleon process 
corresponding to x. The following lemma relates the total amount of ink in this process to 
the near- uniformity of xl conditionally on z^. 

Lemma 8.1 Given x = (z,2;) G {y)k, ^ei x^ = {zf,xi) where, conditionally on z^, xi is 
uniform over V\0{zl). Then: 



dTv(/:[xf]),/:[if])<E 



ink? 



\V\-k + l 



Fil 



where Fill is the event defined in Lemma 7.1 (see also Remark 1.2). 
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Proof: We have seen that zj and zf have the same distribution (cf. the proof of Proposi- 
tion [5]2]). We deduce that: 



,c 



Vt > 0, Vb = (c, b) G {V)k, P {zi = c, ii = b) = ,yY_j^^\ = P ({zf = c} n Fill) 
where the last equahty foUows from Lemma |7.2[ On the other hand, (jS.l.ip imphes: 



Vt > 0, Vb = (c, b) G {V)k, P (xf = b) > E 
We deduce that: 



{zf=c}nFii 



|inkf(5) 



.1.2) 



Vt > 0, Vb = (c, b) G {V)k, (P (zf = c,xl = b) -r {zj = c, xf = b))+ < 



< E 



I{zf=c}nFiii(l-inkf(6)) 



E I|2.c=f.|nFiii(l ~ '"'^f (^)) since the integrand is > 0. 



We now combine this with formula (|2.2.3p for (iTv("i ")• 

dTv(/:[xf],/:[x[]) < Yl IE[l^,c=,}nRii(l-inkf(6)y 

b=(c,fe)G{V)fe 



= E ^ 

ce(v)fc_i 

(sum over b + (l6X2]l ) = J2 ^ 

ce(v)fc_i 



{zC=c}nFiii E (^ ~ '"^^t{b)) 

bev\o{c) 



I{zf=c}nFiii(l^|-^ + l-inkf 



(sum over c) = E []lFiii(|l^| — A; + 1 — ink: 
(apply Lemma Ell = ! 



V -k + l 



D 



8.2 Distance to the stationary distribution in terms of ink 

Proof: [Proof of Lemma 16. Ij We will prove the following stronger inequality: 

inkr 



sup dTv(-C[xf],£[yf]) < 2/i; sup E 
x,ye(v)fc w6(y)fe 



1 



\V\-k + l 



Fil 



?.2.11 



which implies the Lemma by convexity. 
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Declare two states u, v G (^)fc to be adjacent (u ~ v) if they differ at precisely one 
coordinate: i.e there exists an i G [k] with u(i) 7^ v{i) and u(r) = v(r) for r G [A;]\{i}. We 
first bound d^y {Jl[}i.j-], Jl[yj-]) for adjacent x ~ y. 

One can assume without loss of generality that x and y differ precisely at the A:-th 
coordinate. Using the notation from Section 15.41 we write x = (z, x) and y = (z, y) for 
z G {V)k-i and x G y\0(z). Defining Xj = {zj.,xl) as in Section [8.11 and y/^ similarly, we 
see that C[x(] = C[y(] for all t > 0. We deduce that: 

dTy{C[^(],C[yi]) < dT:y{C[^i],C[^i]) + dTY{C[yi],C[yi]) 

+d^y{C[5ci],C[yi]) 
(3rd. term= 0) = dTv(/:[x[], £[x[]) + dTviC[y't]X[ft]) 

inkf 



(use Lemma |8. II) < E 



< 2 sup E 

we(y)fc 



y| -A: + l 
1 



Fill 



inkr 



+ E 



inkf 



y -A: + l 



Fill 



\V\-k + l 



Fil 



^.2.2) 



Now consider x,y G {V)k arbitrary. One can find a sequence {x[z]}[^q C {V)k with 
r <2k and: 

x[0] = X ~ x[l] ~ x[2] '- x[r] = y. 

The triangle inequality gives: 

r 
dTy{C[^i],C[y(]) = d^y{£[^[0]i],£[^[r]i]) < J2dTy{C[^[i - 1][], £[x[^]f ]). 

Applying (J8.2.2J) to each adjacent pair x.[i — l],x[i] gives (I8.2.ip . □ 



9 Depinkings are fast 

The results in this section lead to the key Lemma 16.21 We first show that, in the first 
two phases of the chameleon process ~ a constant color and a color-changing phase -, the 
number of red particles decreases in expectation by a constant factor. 



Lemma 9.1 (Proven in Section 19. Ih Consider a modified chameleon process where one 
drops condition 2. for a pinkening step (cf. Box \5.1\) . Assume also that k < \V\/2, \V\ > 300 
and that the initial state {z, R,P,W) G Ck{V) with \P\ < \R\ < \W\. If the phase length 
parameter T satisfies: 

T>20T|p(2,G)(l/4), 

then: 



E 



1 7?' 



c 

2T_ 



< {l-c)\R\ 
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where c = 1/1000 > 0. 

With this, we will show that the first depinking time has an exponential moment. 

Lemma 9.2 (Proven in Section 19. 2}) Consider a chameleon process (without the modi- 
fication in the previous Lemma) with \V\ > 300 and k < \V\/2, started from an initial state 
(z, R, P, W) G Ck{V) with \P\ = 0. There exists a universal constant K > such that if the 
phase length parameter T satisfies T > 20T|p(2_G')(l/4), the first depinking time Di of this 
process satisfies: 

r ^1 " 

K eKT < e. 
In Section [9.31 we deduce Lemma 16.21 from Lemma 



9.1 Loss of red particles in the two first phases 

Proof: [Proof of Lemma |9.1| Note that there is no depinking at time t = 0, since there are 
less pink particles than white or red ones in the state (z, R, P, W). Finally, the conditions 
on P, W, R and k imply: 

3\W\ > \R\ + \P\ + \W\ = \V\ - k + I > \V\/2 ^ \W\ > \V\/6. 

The interval (0, T] is a constant-color phase where black, red and white particles are 
simply moved around. Lemma 15.21 shows that the state of the process at time T is given 
by: 

iz^,R^,P^,W^) = iIiz),I{R),I{P),IiW)) 

where I = -^(o,r] = ^T is the map obtained from the modified chameleon construction in 
Section [5l We will need the following properties later on: 

Proposition 9.1 (Proven in Section 19.1.11) For all {a,b) G (1^)2 and S,L C V with 
SxLc{V)2, \L\ > |y|/12, 

P ((a, b) G I{S) X I{S)) < P (a G I{S)) (^ + 2'^^ 

P ((a, 6) G I{S) X I{L)) > ^ (1 - 2-^) > ^. 

Remark 9.1 The intuitive meaning of this is that {Rj^,Wji) are close to uniform in terms 
of correlations of "pairs of particles" at the end of the constant- color phase, and this will 
only hold because T = Q (T|p(2,g))- Morris^ original argument for {'L/L'L) could instead 
rely on good estimates for transition probabilities for single-particle random walks. We note 
that we need the negative correlation property in the proof of this proposition. 
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In the time interval (T, 2T) , each time T < Tm < 2T may or may not be a pinkening 
time, depending on whether pinkening condition 1. is satisfied. We will nevertheless consider 
the maps: 

It = I(^T,t], T <t<2T (cf. defn. in Section EH). (9.1.1) 

We emphasize that It does not correspond directly to the evolution of the chameleon process 
in the time interval (T, 2T). Propositions 13.3 1 and 15.1] implv: 

Proposition 9.2 (Proof omitted) {It}T<t<2T is independent from I, and so are all the 

points of the Poisson process {r„}„ in the interval (T,2T) and all markings en,Cn corre- 
sponding to these points. 

We need a new definition before we proceed. Let a £ V he given. Let (pa be the first 
time of the form Tm with T < Tm ^ 2T for which a G em] if no such time exists, let 
(j)a = +00. If (pa < +00, there exists a vertex b £ V such that the edge Cm just mentioned 
has a = I(j,^_{a) and I,j,^_{h) as endpoints immediately prior to time (pa- We set Fa = hm 
that case, or Fq = * if (pa = +oo. The following simple claim is essential to what follows. 

Claim 1 The number of pinkening steps performed in time interval (T, 2T) is at least the 
number of h £ I{^) such that Fa = b for some a € I{R)- 

Proof: [Proof of the Claim] Let b G I{W). Given the rules for color-changing phases (cf. 
Box 15. ip , the particle at that location will move in the time interval (T, 2T) according to 
It until the first time t = Tm & {T, 2T) such that It_ (b) E em and the other endpoint of 
Cm is white (if such a time exists). Now if a E /(-R) satisfies Fa = b and Tm = (pa, we have 
{/t_(6),a} = em and a must still be red at time {(pa)-, since it was not contained in an 
edge before in this phase. It follows that the particle started from b must become pink by 
time (pa- n 

The claim implies: 

l^2r_ I = l-R| - # of pinkening steps in (T, 2T] (9.1.2) 

< 1^1- E ^^..nn){Fa=b}- (9.1.3) 



h(^I{W) 



The sum in the RHS satisfies: 



b&I(W) aeI{R),beI{W) {a,a'}cIiR),b<£l(W) 
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and we obtain: 
^\\R2tJ-\R\\<- Y. ^{a(^I{R),beI{W),Fa = b) 

{a,b)e{V)2 

+ Yl P{a,a' Gl{R),beIiW),Fa = b,Fa>=b). (9.1.5) 

{a,a'}cV,b£V 

The event {Fa = b} is entirely determined by the points of the marked Poisson process 
and by the coin flips performed in the time interval {T,2T), and therefore is independent 
of/ (cf. Proposition I9.2p . We deduce: 



Y r{aGl{R),bGl{W),Fa = b) 

{a,b)e{V)2 

= Yl ^{aeI{R),beI{W))F{Fa = b) 



(a,6)e(V) 



2 



(use Proposition EM > T^L Y ^ i^a = b) = -^ ^ P (F, / *) . (9.1.6) 

{a.,b)&{V)2 ' ' aeV 

For a given a £ V, F {Fa = *) is the probabihty that there is no T < r^ < 2T with en 3 a. 
Notice that this is at most the probabihty that I2t{0') = a: a cannot move if there is no 
edge Sn B a with T < Tn < 2T. We deduce: 

P (Fa / *) > 1 - IP [hria) = a) = P (a^ ^ a) 

where {a^'}t>o is a reahzation of RW(G) started from a. By the contraction principle and 
Proposition 12. H 

T > 20T|p(2,G)(l/4) > 20Trw(g)(1/4) > Trw(g)(2~'°) 

which implies: 

1 13 

P (a^ / a) > 1 - — - 2^20 > _ gij^(,g |y| > 300. 

We deduce from ()9.1.6p that: 

Y r{a€liR),b€liW),Fa = b)>^-^. (9.1.7) 

(a,6)GF2 
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We now consider the second sum in the RHS of (|9.1.5p . As before, we notice that 
{Fa = b, Fa' = b} is independent of / and therefore: 

Y, F{a,a' eI{R),beI{W),Fa = b,Fa'=b) 

{a,a'}cV,b£V\{a,a'} 

Y, F{a,a' £liR),b£l{W))FiFa = b,Fa'=b) 

{a,a'}cV,beV\{a,a'} 

< Yl r{a,a' €l{R))F{Fa = b,Fa'=b). 

{a,a'}cV,beV\{a,a'} 

We claim that: 
Claim 2 (proven in Section [9. 1.2h For all {a,a',b) £ (y)3, 

¥{Fa = b,Fa' =b)<F{Fa = a' , Fa' = b)+F{Fa' =a,Fa = b). 
Summing up over b above gives at most P {Fa = a') + ¥ (Fa' = a) in the RHS. Therefore, 

J2 F{a,a' Gl{R))F{Fa = b,Fa' = b) 

{a,a'}cV,b&V\{a,a'} 

< Y ^{a,a' eI{R)) {F{Fa = a')+F{Fa' = a)) 

{a,a'}cV 

= Yl ^ (« ^ ^(^) ' «' ^ ^(^)) ^ i^<^ = "') 

{a,a')e(V)2 



(apply Prop. [9lD = (^ + 2'^^ Y ^^ (« ^ ^(^)) ^ (^" = "') 

{Ua'iFa = a'} = {Fa / *}) = (^ + '^''') T.^ i^ ^ ^R)) ^ {Fa + *) 

^' ' ' a&</ 

(P {Fa /*)<!) < (^ + 2-^°) 5^ P (a G /(i?)) 

= I — — ■ + 2^ I |i?| since I = -/^(o,t] is a bijection. 
Plugging this equation and (I9.1.7P into (19.1.50 we obtain: 

E[I^?t|-|^I] < 1^1 (^ + 2"'°- n) ^ -|i?|/30, if|iZ| < |y|/28. (9.1.J 
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If \R\ > 1^1/28, we can still find a subset Rq C R of size \Ro\ = Ll^l/28j . Since: 

XI h3aeI{R):Fa=b}> "^ ^3ael (Ro) ■ Fa=b} , (9.1.9) 

b(^I{W) beI{W) 

we may repeat the reasoning presented from (j9.1.4p onwards, replacing Rhy Rq, to deduce 
that: 

^ l-^ol 



E 
We now note that, since \V\ > 300, 



^2T_I - l-^l 
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iR,i>M_i>m>3w 



30 - 100 - 100 

since \R\ < \V\. We deduce that: 



^?tJ - l^l] < -^ if 1^1 > 1^1/28, 



E 
which gives the Lemma together with (j9.1.8p . □ 



9.1.1 Proof of the required estimates for the / map (Proposition [97T]) 

Proof: [ProofofPropositionEHRecanthatr > 20T|p(2,G)(l/4), therefore T > 2T|p(2,G)(2~^°) 
by Proposition l2.1[ By the contraction principle [l], this also implies that T > Trvv(g)(2~^'^)- 
Recall that / = /(o,ri ^^ in the construction of the modified chameleon process. This 
implies that for any set S, I{S) has the law of EX(|5|,G) started from S. We deduce: 

P((a,6) G/(5) x/(5)) = F{{a,b}cS^) 
(negative correlation. Lemma IHTT]) < P (a G S^) P (6 G S^) 
{C[I] = C[I-^], Proposition [32]) = P(a e 1(5)) P (6^ G S) 

(r>TRw(G)(2-^°)) < P(ae 1(5)) (1^ + 2-10 

As for the other inequality in the proposition, we have: 

P((a,6) G/(5) x/(L)) = P((a^,6^) e5xL) 
(take X = (a, b)) = P (x^ G 5 x L) 



> (1-2" ^ 

> (1-2-8) 
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l(^)2| 

|y|2 ' 



where (*) follows from the symmetry of the transition rates of IP(2, G), the fact that T > 
2T|p(2,G)(2"^°) and Proposition O We note that \L\/\V\ > (1/12) and 1 - 2"^ > 12/13 
to finish the proof. □ 

9.1.2 Proof of claim on Fa (Claim [2]) 

Proof: [of Claim [2] It suffices to show that for {a,a',b) G (^)3, 

P {Fa = b, Fa> =b,(t)a<<Pa')=^ {Fa = b, Fa' =a,^a< M ■ (9-1.10) 

Let Lfe, Rb denote the events appearing in the LHS and RHS of (|9.1.10p (respectively). We 
present a simple measure-preserving mapping $, which acts on 

{'P,{en}n,{Cn}n,{dk}k), 

that maps Lh into Rb and vice- versa. We describe $ in words: all values of dk,T < tj < 2T 
and all corresponding ej and Cj, except for the following modification: if Tm = (pa, we flip 
the value of Cm to c^ = 1 — c^ . 

Let us check that <1> has the desired properties. ^ is clearly measure-preserving, since 
(pa is a stopping time that is independent of the value Cm of the flipped coin. 

Now suppose {It}T<t<2T is defined precisely as {It}T<t<2T, but with Cm flipped. It is 
easy to see that (pa,(pa' retain their values and that the random variable Fa corresponding 
to Fa in the / process satisfies Fa = F^.The two processes coincide for any time T < t < (pa- 
If Lb holds, we have (pa = tj < 2T for some j, and the endpoints of ej are a and /<^^ {Fa) = 

Itf,^_ (by definition of Fa). Since the coin fiips used for /<^^ and /^^ are opposite, we have 

(4(a),/^,(6)) = (1^,(5), 4(a)) 
whereas I(f)^{c) = I^^{c) for all c G y\{a, F^}. 

Under L,, \/t G [0„,2r] : {it{a) , it{b)) = {it{b) Jt{a)) . 

Under Lb, the edge e^ corresponding to r^ = (pa> was of the form e^ = {a',I(j, ,_{b)}. This 
implies e^ = {a',I^ /_(o)} in the It process, and the latter must be in the event Rb- This 
shows that ¥ {Lb) < IP {Rb)- The opposite inequality follows from reversing roles of the two 
processes. □ 

9.2 Estimate for the first depinking time (Lemma 19.21) 

Proof: [of Lemma [9. 2 j As in Lemma [9. II we drop condition 2. for a depinking time, and no- 
tice that this change does not change the value (or the distribution) of Di . The modification 
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to the process also does not affect the end result of Lemma 15.11 that is, the discrete-time 
process starting from {z, R, P,W) with subsequent states {zi, Ri,Pi,Wi) described in that 
Lemma is a time- homogeneous Markov chain, and Di = Di/2T is a stopping time for this 
process. 

Moreover, we assume without loss of generality that \R\ < \W\ , which implies that 
|i?p| < \Wf\ for t < Di. This implies \Wi'\ > |V^|/6 unless there are more pink particles 
than red ones at time t < Di; this follows from the reasoning in the beginning of the proof 
of Lemma 19.11 in Section 19.11 Recall that each pinkening step remores a red particle and 
creates two pink ones. It follows that |-R^2-_| < 2|-R|/3 implies Di < 2kT, and: 

yk£-M+,F(Di>k] < F(\Rk\>2\R\/3, Di>k-l) 



< 



3E 


^k^{D,>k} 




3E 


2\R\ 
E [kk 1 fk 


-1 



\Dj>{k-l)} 



2\R\ 



(9.2.1) 
(9.2.2) 



where Tk-i is the a-field generated by (zj, Ri, Pi, Wi) ioi i < k — 1. 

We now estimate the integrand in (j9.2.'2|) . Lemma 15.11 and its proof implies that 



E 
co! 



Rh I J^i 



k-l 



is the expected number of red particles after a potential depinking, a constant- 



or phase and a color-changing phase for a chameleon process started from 



{ik-i,Rk-i,Pk-i,Wk-i) G Ck{V). 
By Lemma |9. 11 we can ensure that: 



E 



Rk I J^k-1 



< (1 - c)Rk-i if \Wk-i\ > \V\/Q and |A-i| < \Rk-i\. 



As noted before, these conditions are always satisfied in the event {Di > (k — 1)}, because 
there are less pink than red particles. We deduce: 



E 



yk G N+, 



-^*:^{Di>fc} 



\R\ 



E 



E 



Rk I ^k-l 



< 



< 1-c 



{Di>ik-1)} 



E 



\R\ 

Rk-1^0^yk-l} 



\R\ 



(...induction...) < (1 - c)^. 
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This implies: 



■(L»i >2kT) =^[bi>k] < 



3(1 - c) 



-, c= 1/1000 universal. 



From this one can easily show that E [e^^'^^] < e for some universal K. □ 



9.3 Proof of Lemma 16.21 

Proof: Fix x G {V)k- We first prove that: 



E 



e KT 



< e-' , K > from Lemma 



(9.3.1) 



this is the bound we wish to obtain except that we are not conditioning on Fill. 
We proceed as in the previous proof and consider the discrete-time process 

{{zi,Ri,P„Wi)}i>o 

introduced in Lemma l5. 11 henceforth called the hat process. This time we take the initial 
state 

(z, R, P, W) = (z, {x}, 0, y\(0(z) U {x})) 

corresponding to x = {z,x) in the sense of Proposition 15. 2i Also recall the definition 
Di = Di{x.)/2T and note that (19.3. ip is equivalent to: 



E 



eK^ 



< e^, K' = 2K. 



(9.3.2) 



This is valid for j = 1 due to Lemma l9.2i For j > 1, we recall the definition of the a-fields 
J-i, recall that -Dj-i is a stopping time for the hat process (cf. Lemma l5.ip and obtain: 



E 



eK^ 



< E 



e~K^E 



e K^ 



-F, 



D 



3-1 



(9.3.3) 



We will apply the strong Markov property of the hat process (cf. Lemma 15.11 again) to 
bound the conditional expectation in the RHS. The conditional law of Dj — -Dj-i given 
J-'a is the law of the hat process started from state: 

Notice that P^ = P£ ^ 0; in fact, since depinking occurs at time Dj_i, we know that 
\Pd I — niin{|i?^. U^/^^. |}. However, at time Dj-i all pink particles disappear: 
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the hat process evolves as if started from a state with no pink particles, and Dj — -Dj-i 



is the first depinking time for the hat process with this modified initial state. 
from Lemma 19.21 that: 



We deduce 



E 



-^j-i 



Tf 



so that: 



E 



CK' 



< E 






< e almost surely, 



e < e-' by induction. 



(9.3.4) 



This proves (I03]l and (|9Xn) . 

To prove the Lemma, notice that conditioning on Fill simply biases the coin flips dk 
performed at depinking times (cf. Remark l7.3p . This will not change the distribution of Di 
or the conditional distribution of Dj 
we presented above still applies. □ 



Dj~i given the past of the process, so the argument 



10 Final remarks 

Our paper leaves many questions open. Here we present a few problems that seem especially 
interesting: 

• Are there any other interacting particle systems whose mixing parameters can be 
bounded solely in terms of the constituent parts? Non-symmetric exclusion is an 
obvious candidate. Another is the zero-range process. Morris jl8] used the comparison 
principle and a coupling argument on the complete graph to bound the spectral gap 
of this process on a grid. Can one do something less indirect over an arbitrary graph? 

• Can we find a mixing time upper bound of IP(|y|,G) (i.e. as many particles as 
vertices) that is similar to our main Theorem? Inspection of the chameleon process 
shows that it gives the conditional distribution of a particle given the whole past 
trajectory of the other particles. This means, in particular, that it cannot deal with 
A; = |y| particles. 

• Recall the heuristic assumption in the introduction: T^x{k,G) i^) ^ Ci Trvv(g) (^ /^) 
with Ci > universal. Is this actually true? This would be stronger than our main 
Theorem. 

• Combining the previous two items: is it true that rip^yi Q)(e) = Ci Trw(g)(£/I^I)'^ 
Could it even be possible that r|p(|y| Q)(e) < Tf(\N{\v\,G)i^)j i-^- ^^^ interchange process 
mixes at least as fast as independent random walkers? This would give Aldous' (now 
proven) conjecture on the spectral gap as a corollary. 
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A Appendix: mixing bounds for EX(/c, G) via canonical paths 
and log Sobolev inequalities 

We use asymptotic notation below as shorthand; see eg. [2] for precise definitions. Let 
G = {V,E,{we}eeE) be a weighted graph. It seems that the best general bound that was 
previously available (implicitly) for the mixing time of EX (A;, G) comes from the combination 
of three ingredients. 

Mixing time from Log-Sobolev constant. The state space of EX{k,G) has car- 
dinality Cfc') = 2®(l^l) if A; = e {\V\). By the results of [9], if PEX{k,G) is the log-Sobolev 
constant of EX{k,G), then: 



TEX(fc,G)(l/4) = O f±}}l-) foTk = e {\V\ 



\PEX{fc,G), 

Log-Sobolev inequality for the Bernoulli-Laplace model. Consider the complete 
graph Kiyi where each edge has weight 1/|V^|. EX{k, -f^ivi) is the so-called Bernoulli- Laplace 

model with k particles, whose log Sobolev constant is of the order Q ( ln(|y| /A;(|y| — k)) j . 
Notice that this is 9 (1) for A: = G {\V\) 

Comparison argument. Now consider a general weighted graph G = {V, E, {we}e£E)- 
Assume that for each pair (x, y) £ V'^ one has defined a path 'jx^y in G connecting x to y. 
For each such pair, let Ix,y{e) = 1 if e is crossed by jx,y and otherwise, and also let £x,y 
denote the length of 'yx,y Finally, define: 

,^(G) = max V %ZpM. 

(x,y)£V-' 

It is shown in the proof of |8, Theorem 3.1] that this comparison constant for the Dirichlet 
forms of RW(G) can be "lifted" with no loss to EX(A;, G). The comparison principle for the 
log Sobolev constant [9] implies PEX{k,G) = ^ {4'"^{G))- We deduce: 

TEX(fc,G)(l/4) = O (0(G) In \V\) ifk = e {\V\) ■ (A.0.5) 

It can be very hard to find good upper bounds on 0(G) in general, but the general lower 
bound we will present implies that: 



2dist2 



0(G) > ^-, (A.0.6) 

a 
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where dist is the average over all {x, y) G V"^ of the square of the graph-theoretic distance 
between x and y, and d is the average (weighted) degree in G. Indeed, it suffices to see 
that: 






We 



^l^feE^f 



E 



Jx,y)€V^ 






{nse ^^Ix,yie) = ix^y) = 
(use 4,?; > dist(x,y)) > 



E/e£ ^/ 



1 



E 

.{x,y)£V^ 



\y\ 



1 



E ^x,y 
|\/|2 
{x,y)eV^ ' ' 



E 

.(a;,y)ey2 



dist(x,y)' 

\V\' 



2dist' 



d 



We note that this is a lower bound, which we do not know how to achieve in the examples 
in Table [H 



B Appendix: the trajectory of inkj given Fill 

We use the facts proven in Section [7] to derive the technical estimate in Proposition 16.11 in 

the proof of Lemma 11.21 (cf . Section 16. ip . 

Proof: [of Proposition 16. Ij We take the notation in Section [7] for granted, but omit the 

superscript w in this proof. Our first goal will be to show that, conditionally on Fill, 

{inkj}j>o is still a Markov chain. Repeating the steps of the proof of Lemma l7.ll we note 

that: 

E 



P Fill I (ink 



i)i<j 



inkoo I (inkj)i<j 



inki 



\V\-k + l \V\-k + l 

We deduce from Bayes' rule and the Markovian property that: 

P ( f]{UAi = ai} I Fill ) = F ( fjii^ki = Gi} ) aj 



(Fill) inkj 



\i=l 



',1=1 



(Markov property for inkj) 



p{l,ai)p{ai,a2) ■ . .p{aj^i,aj) aj 
q{l,ai)...q{aj^i,aj) 



where 



bp{a, b) 
q{a, b) = it a 7^ 0. 
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Notice that, since ink^ does not visit in the event Fill, we do not need to define q{a, b) for 
a = 0. We have shown: 

Proposition B.l Conditionally on Fill, the trajectory of {inkj}j>o is that of a Markov 
chain in {1, . . . , \V\ — k + 1}, with transition rates q{a, b) and started from inkg = 1. 

For the remainder of the proof, we wih use this Proposition to bound 1 — ink^/(|y|— /c+1). 
ActuaUy, another quantity is easier to bound. Set Ig = ink^/(|y| — A; + 1) and: 



Zi 



Y^minjl - Ii,Ii} 
h 



Notice that conditionally on Fill, /^ > always, hence Zi is a.s. well defined for all L 
Moreover, one can check that 1 — I(^ < Z^ always. Therefore the Lemma will follow from 
the estimate: 

E"^"' [Zi] < (71/72)Vl^|-fc + l, 

where E*^'" [•] corresponds to an expectation with respect to the conditional distribution 
given Fill. Since Zq = \/\V\ — /c + 1, the above estimate follows directly from the following 
claim. 



Claim 3 



Fill 



V£ G N, E""'" [Zi\ < (71/72) E""'" [Z^_i] . 



Therefore, proving this Claim will finish the proof. 

To prove the claim, we first note that for all i Zi is a function of inkj, and Z^-i = 
Zf = 0. We deduce: 



E""'" [Zi] = E 



Fill 



E 



Fill 



^e-1 



ink£_i 



Zi-iI{z,_,=io} 



(B.0.7) 



We now bound the conditional expectation in the RHS. We may assume that ink£_i = r 
with < r < \V\ — k + 1 (otherwise Z^^i = 0). Thus we wish to bound: 



E^ 



Zf, 



z 



ink^_i = r 



i-i 



,l<r<\V\-k. 



If we note that: 

Zi _ ^min{l - liJi} h-i _ A/min{l - hJi} ink<._i 



Zi 



-1 Y^minjl - h^iJi^i} h i/min{l - /£_i,/<?_i} ink^ 
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and define /(a) = y^minla, |y| — fc + 1 — a}, we see that: 



E' 



Zi 



Z^ 



ink 



t-\ 



i-i 



E 



Fill 



/(ink£) inkf_i .-— - 
— :^^:^ X ^^^ — I ink£_i = r 



/(ink^_i) ink£ 



(use Proposition IB. ip = N^g(r, s] 



/(g) ^ ^ 

/(r) s 



(use formula for q{-,-')) 



fir) 



where p(-, ••) are the transition rates of the unconditional {inkj}j>o process. Using the 
formula for these we obtain: 



E 



Fill 



Z, 



ink 



i-i 



£-1 



1 / /(r + A(r)) + /(r-A(r)) 

2 V f{r) 



(B.O.J 



Recall the formula for A(r), cf. Proposition 17.31 

rminjr, \V\ — k + \ — r\ 

A(r) = 

^ ' 3 

We now split the analysis of the RHS of this in two cases. 

Case 1: 1 < r < (|y| - A; + l)/2. In this ca se /(r) = ^/r and A(r) = [r/3] > r/3. We 
use the upper bound f(r it A(r)) < ^Jr it A(r) to obtain: 

Zp 



w 



ink/> 



1 



A(r) 



1 + 



A(r) 



(B.0.9) 



Recall the bound 'Vl - 2; + \/rTx < 2(1 - x^/S)", valid for ah < a; < 1; this can 
be checked by squaring both sides of the inequality. In our case, we apply this with x = 
A(r)/r > 1/3 and deduce: 



E^ 



Zo 



Z> 



ink£_i = r 



i-i 



,_ir^v<!i. 



(B.O.IO) 



Case 2: {\V\-k + l)/2 <r < \V\-k. In this case (IB. 0.911 holds with r' = \V\-k + l-r 
replacing r. Similar calculations imply that the conditional expectation is also < 71/72 in 
this case. 

Thus we see that in both cases: 



E^ 



'■^i-i 



ink^_i = r 



< 



71 

72' 



Plugging this into ()B.0.7p gives: 

71 

72 



E^'" [Ze] < ^ E^'" [^£-iI{z,_,^o}] = ^ E"" [Zi-i] 



72 



which finishes the proof. □ 
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